Dynamic scaling of fronts in the quantum XX chain 
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The dynamics of the transverse magnetization in the zero-temperature XX chain is studied with 
emphasis on fronts emerging from steplike initial magnetization profiles. The fronts move with 
fixed velocity and display a staircase like internal structure whose dynamic scaling is explored both 
analytically and numerically. The front region is found to spread with time sub-diffusively with the 
height and the width of the staircase steps scaling as t~^^^ and t^^^ , respectively. The areas under 
the steps are independent of time, thus the magnetization relaxes in quantized "steps" of spin-flips. 
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Fronts often emerge in relaxation processes. Sim- 
ple examples are the domain walls generated in phase- 
separation dynamics 0] but there is a long list of fronts 
(also called shocks, active zones, reaction zones, etc.) 
resulting from unstable dynamics in various physical, 
chemical and biological systems 0, Q- Furthermore, 
fronts can also emerge due to the initial spatial sepa- 
ration of stable and unstable states 4:]. 

Apart from being the instruments of relaxation, the 
importance of fronts is also due to their "catalytic" na- 
ture. Namely, structures are built in them 5] and pat- 
terns often emerge in the wake of moving fronts 6] . Thus 
it is not surprising that much effort has been devoted to 
the description of their spatio-temporal structure 

In contrast to fronts in classical systems, not much 
is known about quantum fronts. The examples we are 
aware of are restricted to quantum spin chains: fronts 
or shocks have been seen in the XX model 0, 0j 
in the transverse Ising chain [lol as well as in the 
Heisenberg model (1^ [13 ■ The detailed structures of the 
front regions, however, have not been elucidated even in 
these cases. 

Having in mind the importance of the front dynam- 
ics, we set out to investigate the fronts emerging in the 
quantum XX chain. The transverse magnetization in 
this model is conserved and its relaxation is known to 
be governed by fronts provided an initial state with a 
spatial separation of distinct magnetizations (±mo) is 
prepared [g. The problem with such initial condition is 
exactly solvable and thus we can follow the evolution of 
the spatio-temporal structure of the front in detail. 

The main result of our calculations is that, in addition 
to the known global scaling (finite front velocity and 
well defined magnetization profile in the t oo limit), 
there exists a dynamic scaling regime in the front, and the 
associated scaling function exposes a staircase structure 
in the magnetization profile. An important feature of the 
steps in the staircase is that while their height decreases 
as t~^/^ and their width increases as i^/'^^ the areas under 
the steps are constants and, furthermore, the constants 
are the same for all the steps. The value of the constant 



(twice the magnetic moment of a spin) indicates that a 
step carries a reversed spin with respect to the aligned 
initial state. These reversed spins move with the front 
velocity, keep their identity with respect to other reversed 
spins, and their spatial spread is sub-diffusive. Thus the 
magnetization relaxes in well defined quantized (in space 
and in time) "steps". 

Before turning to the calculations, we note that our re- 
sults suggest that the quantum fronts and, in particular, 
the steps carrying a unit of spin-flips, may be envisioned 
as ingredients in controlled transport of bitwise informa- 
tion in magnetic nanostructures. 

The system we investigate is the XX chain defined by 
the Hamiltonian 



H = - 



N-l 

E 



(1) 



where the spins S" {a = x,y, z) are 1/2 times the Pauli 
matrices situated at the sites [j — 0, ±1, ±(A^ — 1), A^] 
of a one-dimensional lattice. Free boundary conditions 
are used, and the thermodynamic limit A'^ — > cxd is as- 
sumed throughout the paper. The physical quantities are 
measured in their natural units: energy in Jti^ where J is 
the nearest-neighbor coupling, magnetization in h, length 
in units of the lattice constant a, and time in 1/Jh. 

Our principal aim is to investigate the time evolution 
of the (globally conserved) transverse magnetization 



m(n,t) = (^|5^(t)|^) 



(2) 



emerging from an initial state | ip) with a steplike mag- 
netization profile 
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Since the dynamics following from H can be described 
[1^ lisj in terms of local fermionic operators (c„,c]^) 
whose Fourier transforms diagonalize H, and since 
can be expressed through the local fermionic operators 
as — cl^Cn — 1/2, the evaluation of ((^| S^{t) \ip) is a 
relatively simple exercise. The calculations have been 
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carried out in [3 with the result for n> 1 given through 
the Bessel functions of the first kind 



m{n,t) = -ijoW - ^ 
1=1 



J? it) 



(4) 



while the expression for n < is obtained from symme- 
try considerations m{n,t) = —m{—n + \,t). The global 
scaling of m(n, t) emerges in the n — > +00, t — > +00 and 
n/t = finite limit where m(n, t) can be written in a scal- 
ing form m(n, i) — > $ ( j) , and the scaling function $ (v) 
is given Q by 

{— TT^^ arcsin(u) for < w < 1, 
(5) 
-1/2 fort;>l. 

The Too = 0.5 curve in FigQ shows together with 
the shape of m(n, t) plotted against v — n/t at finite time 
{t = 200). As one can see, the magnetization displays a 
staircase structure near the edge of the front {n/t « 1). 
Our main concern will be the scaling properties (both in 
space and time) of this staircase structure. 




FIG. 1: Magnetization profiles emerging from steplike initial 
conditions [eqs. and ©]. The large-time scaling limit is 
shown by dashed lines [eq. (|SJ for mo = 0.5]. 

For possible applications, it may be important that the 
front properties were tunable. Possibility for tuning can 
be seen by solving the problem for more general steplike 
initial states 



(V^mo I S^iiO) I ipmo ) = 



Too 

-Too 



for 
for 



- < n < 
1 < n < iV . 



(6) 



where | (pmo ) is constructed by joining two half chains 
which are the ground states of the XX model at magneti- 
zations ±TOo. The expression to be analyzed for this case 
can be taken from Ref.0, eq.(12). It is more involved 
but its numerical evaluation does not pose difficulties. 
We shall not describe the numerical work but results for 
Too 7^ will be discussed and displayed (Figsnand|31). 



The present work on the dynamic scaling of the mag- 
netization profile around the edge of the front (n « t) 
originates from our numerical studies of the deviation of 
the magnetization 6m{n,t) — m{n,t) — m(t,t) from its 
front value m{t,t). As one can observe on Fig|21 a well 
defined scaling function emerges in the limit < — *■ 00 pro- 
vided the magnetization and the region around n « t are 
scaled by i^/^ and t~^/^, respectively. 




z = (n-t) / 



FIG. 2: Scaling in the front region. The magnetization mea- 
sured form m,(t, t) is magnified as t^^^ ■ (5m(n, t) while the dis- 
tance from the front position is scaled as z = {n — t)/t^^^. The 
analytically derived scaling limit [ea. lllll ] is indistinguishable 
from the t — 10^ curve in the 2-range plotted. 

An analytic understanding of the scaling seen in Fig|21 
can be developed by first using Q to write Sm{n,t) as 



Sm(n, t) = < 



i=t 



for n < t, 
for n ~ t, 
for n > t. 



(7) 



Next, one can observe that, in what we call the scal- 
ing regime {t — + 00, n — ^ 00, and \t — 7i| « t^^"^), 
the sums above contain Bessel functions whose index 
(I) and argument (t) differ from each other at most by 
\l — t\ ~ ©(t^/"^). This suggests the use of the following 
asymptotic expansion of Bessel function [l^ : 

J„(n + zn^l"^) « 2i/3„-i/3Ai(_2i/3^) + 0(^-1) (8) 

where Ai(z) is the Airy function. Indeed, in the scaling 
regime, the terms in the sums in eq.Q can be written as 



22/3 
7273- 



Ai^ 2^/3 



,l-t 

7W 



1 



f2/3 \ ^1/3 



l-t 



(9) 



where F{y) = 2"^/^ ki^{2^/^y). Note that the last expres- 
sion in (PJ is obtained by replacing Z by i in the factors 
?^/3 and P/^^ an approximation that is valid in the scaling 
regime. 
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The last step in deriving Sm{n,t) is the calculation 
of the sums in lO using eq.©. In the scahng regime, 
the sums can be replaced by integrals and, introducing 
the variable y = {I — t)/t^^^, both sums yield the same 
expression. As a consequence, we obtain a single scaling 
form for n > t as well as ior n < t 



5m[n,t)^^G 



(10) 



where the scaling function G{z) is explicitely given by 
the following integral 



dyF{y)^-2^/' dyA\'{2'/'y). (11) 



G(z) 



The scaling form IjlUII and the scaling function (|ll|l are 
our central result. The staircase shape of the front fol- 
lows from (|llll while the scaling form (|10|1 shows that 
the width of the steps increases a.s w ^ t^^^ while their 
heights (measured from the ^ ^ oo level) decrease as 

An important feature of the above scaling is that since 
w X h ^ t'^, the areas under the steps are constants. 
Thus a step carries a given amount of magnetic moment 
Ids = Ws X hs where subscript s denotes the s"^ step 
counted from the edge of the front. In principle, fis could 
depend on s but their numerical evaluation for the steps 
near the edge of the front (1 < s < 8) (see Table Hj) 
strongly suggests that — ^ (the deviations from 1 can 
be attributed to the finiteness of the sums at finite t). 



Ai(2i/3z) 
s'^ zero is given 



equation. The large s asymptotic of the 
13 by \zs\ « (37rs)2/3/2 and thus 



K = [3s/(^2^)]i/3 



(12) 



We need to find now the width of the s*'^ step. Defin- 
ing the borders of the steps as the consecutive inflection 
points on the staircase, one can see from (lllf) that the 
width of the step is given by the consecutive zeros of 
Ai'(2^/3zs) = 0. The large s asymptotic of these zeros are 
known again |Q| and we find \zs+i — Zs\ ~ [7r^/(3s)]^'''^. 
We should remember now that hs was calculated for the 
original (unsealed) magnetization, and the width of the 
step must also be obtained in unsealed spatial coordi- 
nates. Consequently, we should scale |zs+i ~ by t^/"^ 
and thus 



^ t^/^\Ss+i - Ss\ 



(13) 



Comparing eas. ((T^ and ((T^ we see that hs • Ws = 1 
thus arriving at an important property of the staircase 
structure, namely, the steps of the staircase carry a unit 
of magnetization. 

In order to see the relevance of the above result and to 
develop a picture about it, let us consider the magnetiza- 
tion flux in the front region. The total transverse magne- 
tization = J2n is conserved and so one can define 
the local magnetization flux, = SyS^_^_i — S^S^_^_i, 
which is related to the local magnetization through the 
continuity equation. Thus, not surprisingly, one can de- 
rive a simple expression for the time evolution of the 
expectation value of as well 



step n°:s 


A*s(t=103) 


/is (1=10-*) 


Ms(t=10^) 


1 


1.0183 


1.0420 


1.0379 


2 


0.9707 


1.0154 


1.0094 


3 


1.0371 


1.0080 


0.9934 


4 


0.9418 


0.9710 


1.0124 


5 


0.9022 


0.9970 


0.9986 


6 


0.9588 


1.0070 


0.9860 


7 


1.0123 


0.9434 


1.0128 


8 


0.9286 


0.9866 


1.0037 



TABLE I: Magnetic moments carried by the steps (area below 
the steps in the magnetization profile) at finite times {t — 
10^, 10", 10^). The borders of the steps were defined by the 
infiection points in the magnetization curve. 

The = 1 result can be derived for large-order (s ^ 
1) steps. Indeed, as can be seen from Fig^ the height 
of the step can be estimated from the global scaling 
function TO(n, f) = — 7r~^ arcsin (1 — z/i^/"^) « —1/2 

-1/3 



and thus dm = h^ — 



where Zs is the scaling variable at the s*** step. Its value 
is found by finding the s"^ zero of G'{z) which, in turn 
[see ea. (|ll|l ]. is obtained from the s"^ solution of the 



j{n,t) ^ {ip\JUt)\v) = Y.Mt)Ji+iit). 



(14) 



l—n 



The analysis of Sj{n,t) = j{n,t) — j{t,t) parallels that 
of 5m(n, t) and yields the same scaling structure in the 
front 



n — t 



5min,t). (15) 



Adding G{oo)/t^^^ to both sides of the above equations, 
and taking into account that j{n — > oo,t) — and 
■m{n — > oo,t) = —1/2 and, furthermore, remembering 
that the front moves with velocity w = 1, we can write 
the above equations in a form that is easy to interpret 



m{n,t) - (--) 



(16) 



Since the sum of [m{n,t) — (—1/2)] for an interval gives 
the number of up-spins in that interval, the meaning of 
the unit area under under the steps is as follows. The 
steps represent spatial intervals in which a single up-spin 
is spread in the sea of down-spins. These up-spins be- 
have like particles, they move with velocity v = I and 
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provide the magnetization flux. Furthermore, they are 
locaUzed in the sense that their spatial spread is sub- 
diffusive (~ t^^^) instead of the usual quantum mechan- 
ical spreading (~ t^/^). The picture thus emerging is 
somewhat reminiscent of the hard-core bosons descrip- 
tion of the XX model near its critical external field ■ 
From the point of view of possible applications, it is 
important that the staircase structure can be observed 
for other initial conditions as well. Using initial states 
\ipmo) with ±TOo steplike profile as described in eq.®, 
we found numerically that the staircase emerges again 
(see Figl^))- Furthermore, the steps of the staircase were 
found to have the same size independently of values of 
mo (see inset in Fig. O. What is varied with itiq is the 
number of steps N{t,mQ) in the staircase. 
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FIG. 3: Magnetization profiles a,t t — 10"^ for different initial 
values mo [eq.JSJ]. The inset with the curves shifted vertically 
down (m = m — 0.5 + mo) demonstrates that the steps in the 
front region are independent of the value of mo. 

We can estimate N{t,mo) by assuming (on the basis 
of numerical solutions) that scaling extends to the whole 
front region, up to the point where m reaches the steady 
state value of m = 0. Then N{t, mo) is found by locating 
the A^"* step which has a width wn such that a single 
spin flip changes the magnetization from —mo to 0. Since 
Too ■ wn is the difference between the number of up- and 
down-spins in the N*^ step, the above consideration gives 
the condition — toq • wn -1-1 = 0. For large N, ea. H13|) 
gives the width wjv, and we obtain 

N{t,mo) fun^mlt/S. (17) 

It should be noted that ea. l|13|) is an excellent approxi- 
mation also for step-numbers of the order of unity, and 
so ea. p7|) gives N{t,mo) for small N, as well. Indeed, 
e.g. ea. lfT7|) yields iV « 3 for t = 10^ and too = 0.1 while 
the corresponding curve on FigOl displays N fa i steps. 

An important feature of ea. lfTTjl is the strong depen- 
dence of the step number on the initial magnetization 
(A'^ ~ m^). This gives a sensitive control over how many 



steps arrive at a given point at a given time. Thus, in 
principle, one can imagine applications with the steps 
transferring bits of information in magnetic nanostruc- 
tures. The aim of the present study, however, was only 
to draw attention to the remarkable features of quantum 
fronts. Applications would require answers to a number 
of non-trivial questions. First, can fronts with proper- 
ties found in the XX model be observed in other spin 
chains? In particular, are they present in non-integrable 
systems? Second, how does the step structure change 
at small but finite temperatures? Finally, are the front 
properties robust enough to survive small perturbations 
arising from impurities? These are difficult problems but 
the available analytical and numerical techniques may be 
sufficient to tackle them. 

In summary, our calculations exposed a dynamic scal- 
ing regime in the motion of fronts in the quantum XX 
model. The fronts display a staircase magnetization pro- 
file which has a simple interpretation in terms of sin- 
gle spin-flips spread over the spatial extent of the steps. 
Whether these fronts have a direct application is an open 
question but their properties are intriguing enough to 
look for similar structures in more complex spin-chains. 
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